arXiv: 1503.05716vl [quant-ph] 19 Mar 2015 


Equivalence of matrix product ensembles of trajectories in open quantum systems 

Jukka Kiukas and Madalin Gu^a 

School of Mathematical Sciences, University of Nottingham, Nottingham, NG7 2RD, UK 
Igor Lesanovsky and Juan P. Garrahan 

School of Physics and Astronomy, University of Nottingham, Nottingham, NG7 2RD, UK 

(Dated: March 20, 2015) 

The equivalence of thermodynamic ensembles is at the heart of statistical mechanics and central 
to our understanding of eqnilibrium states of matter. Recently, it has been shown that there is a 
formal connection between the dynamics of open quantum systems and the statistical mechanics in 
an extra dimension. This is established throngh the fact that an open system dynamics generates a 
Matrix Prodnct state (MRS) which encodes the set of all possible quantum jump trajectories and 
permits the construction of generating functions in the spirit of thermodynamic partition functions. 

In the case of continuous-time Markovian evolution, such as that generated by a Lindblad master 
equation, the corresponding MPS is a so-called eontinuous MPS which encodes the set of continuous 
measurement records terminated at some fixed total observation time. Here we show that if one 
instead terminates trajectories after a fixed total number of quantum jumps, e.g. emission events 
into the environment, the associated MPS is discrete. This establishes an interesting analogy: 

The continuous and discrete MPS correspond to different ensembles of quantum trajectories, one 
characterised by total time the other by total number of quantum jumps. Hence they give rise to 
quantum versions of different thermodynamic ensembles, akin to “grand-canonical” and “isobaric”, 
but for trajectories. Here we prove that these trajectory ensembles are equivalent in a suitable limit 
of long time or large number of jumps. This is in direct analogy to equilibrium statistical mechanics 
where equivalence between ensembles is only strictly established in the thermodynamic limit. An 
intrinsic quantum feature is that the equivalence holds only for all observables that commute with 
the number of quantum jumps. 

PACS numbers: 05.30.Rt,05.30.-d,64.70.P- 


I. INTRODUCTION 

There is currently much interest in achieving a funda¬ 
mental understanding of the non-equilibrium dynamics 
of quantum systems in general m, and of open quan¬ 
tum systems in particular HHt]. When the dynamics of 
an open quantum system is, to a good approximation, 
Markovian (i.e. memoryless), its evolution is given by 
a master equation of the Lindblad type [S]. Such dy¬ 
namics can be “unravelled” in terms of stochastic quan¬ 
tum trajectories Each quantum trajectory rep¬ 

resents a different observation record (as measurements 
on the environment and consequently on the system, via 
the system-environment interaction). Averaging over the 
whole “ensemble” of dynamical trajectories at any point 
in time gives the average system evolution described by 
the master equation [3 EQ] . 

The concept of trajectory ensembles is important as 
it allows to connect [nns] open quantum dynamics to 
ideas of ensemble statistical mechanics [HI [15] . Consider 
dynamics between zero initial time, up to some time t. 
The dynamics generated by a Lindblad master opera¬ 
tor is that of a continuous quantum Markov chain, one 
of quantum jumps occurring stochastically with dissipa¬ 
tive, yet deterministic, evolution between the jumps. The 
ensemble of trajectories is the set of all possible trajecto¬ 
ries generated by this dynamics, each with an associated 
probability (or amplitude) of occurring and classified by 
one (or more) parameters that are held fixed, in this case 


the total time extent t of the trajectories. In a thermo¬ 
dynamic analogy (and for simplicity we use the language 
of classical systems) if we think of a (canonical) ensemble 
in a system of fixed volume at say infinite temperature, 
the ensemble is the set of all possible configurations of 
the system compatible with that volume, all occurring 
with equal probability (due to the infinite temperature). 

In thermodynamics a central role is played by order 
parameters and their conjugate fields [H El]- Order 
parameters are system extensive observables that quan¬ 
tify properties of the system at the given conditions (e.g. 
density, magnetisation, etc.) which couple to their cor¬ 
responding intensive fields (chemical potential, magnetic 
field, etc.). In the case of quantum trajectories, order pa¬ 
rameters are time-integrated quantities [T21E3], the num¬ 
ber of quantum jumps K over a whole trajectory being 
one typical example. Since the trajectories are stochastic, 
a quantity such as K will be random, and its distribution 
will encode the statistical properties of the dynamics. In 
particular, interesting properties may be uncovered by 
considering atypical values of the order parameter, since 
often significant fluctuation behaviour is encoded in the 
tails of their distributions [T2lE3|. 

This leads then to the idea of “biased” or “tilted” en¬ 
sembles, i.e. ensembles conditioned on certain values of 
the order parameter [I1E3ES]. In the standard ther¬ 
modynamic case, if for example the order parameter is 
the energy E, the ensemble of configurations with some 
fixed E all with the same probability is termed micro- 
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canonical, while the ensemble of all possible configura¬ 
tions with probability proportional to where T is 

a fixed (scaled) temperature, is termed canonical. The 
same can be done in the dynamical case by defining tra¬ 
jectory ensembles conditioned on, say, a certain fixed 
number of emissions K (analogous to the microcanon- 
ical ensemble), or where K is not fixed but trajectories 
are biased by a factor where s is a so-called count¬ 

ing held. This “thermodynamics of trajectory” approach 
[laila], sometimes called s-ensemble method [IlHIS], al¬ 
lows to uncover key properties of the dynamics—often by 
means of large-deviation techniques m that become ap¬ 
plicable when time is large—for example by revealing the 
existence of dynamical phase transitions [laiisiiiiHis]. 

A fundamental property of thermodynamic ensembles 
is that of equivalence: for a given system the microcanon- 
ical and canonical ensemble are equivalent in the limit 
of large system size, as long as the temperature of the 
latter is such that the canonical average of the energy 
coincides with the microcanonical hxed value. The same 
occurs with other ensembles, such as between those of 
fixed volume and fixed pressure. The simplest manifes¬ 
tation of this equivalence is in the Legendre transform 
that connects the corresponding free-energies, i.e., a cor¬ 
respondence at the large-deviation level. For dynami¬ 
cal trajectory ensembles, it has been shown Ennu that 
for classical stochastic dynamics the microcanonical (in 
the sense of fixed time-integrated observable) and canon¬ 
ical (fixed conjugate counting field) ensembles are indeed 
equivalent, as in the thermodynamics case. Furthermore, 
in Ref. |22| . a new ensemble was introduced, one of tra¬ 
jectories have a fixed value of the time-integrated observ¬ 
able, such as the number of quantum jumps, but time is 
allowed to fluctuate (termed x-ensemble). This is like 
the dynamical version of a fixed pressure (and thus fluc¬ 
tuating volume) ensemble, with x playing the role of a 
pressure (on time). In |22j the correspondence between 
the s-ensemble and x-ensemble was proved at the level of 
large-deviation rate functions, both for the classical and 
open quantum cases. 

In this paper we unify and extend these ideas about 
thermodynamics of trajectories and dynamical ensemble 
equivalence by using the technology of matrix product 
states (MPS). We exploit the fact that the set of quantum 
trajectories of an open quantum system can be encoded 
in an MPS, and show how to do this for ensembles of tra¬ 
jectories with fixed number of jumps K. By exploiting 
the MPS construct we prove in general the equivalence of 
microcanonical and canonical dynamical ensembles, thus 
generalising the classical results of HOI mi mm] to the 
quantum realm. We also prove an equivalence result be¬ 
tween the s-ensemble (fixed time) and the x-ensemble 
(fixed number of quantum jumps); in this case we con¬ 
sider the overlap between the pure (system together with 
output) states of the two corresponding canonical en¬ 
sembles, and show that this overlap has subexpoential 
decay in the limit of large time. However, since in the 
time ensemble, the system has a nontrivial evolution af¬ 


ter the last jump, we need to condition onto trajecto¬ 
ries where the last jump happens at the end of the time 
interval. In an alternative approach, we compare the 
(mixed) states of the outputs of the two canonical en¬ 
sembles, rather than the full system and output states; 
we find that the reduced states which describe the output 
up to a fixed time tq converge in the limit of large time 
and respectively counts, and provide explicit expressions 
of the limits. The equivalence is manifest here in that 
the x-ensemble limit is equal to the diagonal part of the 
s-ensemble limit, with respect to the decomposition of 
the output Fock space into “layers” of fixed number of 
counts subspaces. However, the reduced state of the x- 
ensemble does not have any coherence between different 
layers, which may be present however in the s-ensemble. 
We illustrate this purely quantum feature by considering 
the dependence of the two outputs on specific dynami¬ 
cal parameters; in particular the s-ensemble is shown to 
be sensitive to phase transformations generated by the 
counting operator, while the x-ensemble is left invariant 
due to the above mentioned lack of coherence. 

The paper is organised as follows. In Sect. H we set 
up the formalism describing open dynamics in terms of 
MPS. In order to maintain as much generality as possi¬ 
ble we discuss open systems initially in terms of continual 
quantum measurements. This allows to capture both the 
standard situation of open dynamics observed for a fixed 
total amount of time, and also the case where observa¬ 
tions are for a fixed number to quantum jump. While 
the former is known to be encoded in so-called continu¬ 
ous MPS for the system and output, we show that the 
later is encoded discrete MPS states, and describe the 
connection between the two. In Sect. HI we describe 
the thermodynamics of trajectories approach, and discuss 
both the s-ensemble and x-ensemble methods in the con¬ 
text of MPS for the system-output. Section IV contains 
our key results. There we prove the equivalence of en¬ 
sembles. We provide the proof in three stages. The first 
one is the correspondence of potentials, or large-deviation 
rate functions, which duly follows from relations between 
the tilted superoperators that define the ensembles. The 
second stage involves proving the actual equivalence of 
measures over quantum trajectories, in analogy to what 
is done in the classical case. The final stage is the proof 
of full operational equivalence between the ensembles, on 
the level of diagonal blocks with respect to Fock layers. 
In Sect. V we discuss the specific case of quantum renewal 
processes, where the fixed number of jumps ensemble be¬ 
comes particularly simple and transparent, and illustrate 
it with an example. Finally in Sect. VI we provide our 
conclusions and outlook. 


II. CONTINUAL QUANTUM MEASUREMENTS 
AND MATRIX PRODUCT STATES 

We consider the scenario where a quantum system with 
Hilbert space H is continuously monitored by performing 
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measurements on it (typically realised via coupling the 
system with an external environment.) In the counting 
case which we describe below, the measurement record 
consists of clicks of different types, occurring at random 
times, according to a certain probability distribution 
(see also [MIET] for a viewpoint closer to ours). 

Suppose that the system starts in an initial state 
Ip (z %■ For a random waiting time t^ the system 
evolves continuously, so that its unnormalised state at 
time t < tyj is exp{—iHegt)'ip, where Hgg is a non- 
selfadjoint effective Hamiltonian. The probability den¬ 
sity of the waiting time is given by the loss of wavefunc- 
tion normalisation p^(tu,) := —At 
the time of detection, a click with label i = 
is recorded, and the system’s conditional state is up¬ 
dated by applying a jump operator Li. The full mea¬ 
surement is therefore described by the positive operator 
valued measure (POVM) [t^H) ^ J[tuj, i], where 

J\tyj,i\ = is the jump operator effecting the 

total change ft ^ J[tyj,i]'p. These must satisfy 


Nl 

Ptp(.tw) — 'P ' li*]'0ll ) 
i=l 


pOO 

I dt^p^(t^) = 1, (1) 

Jo 


the first condition stating that a detection event happen¬ 
ing at waiting time t^, necessarily involves one of the out¬ 
comes i = 1 ,. .., Nl being realised, and the second fixing 
the total normalisation of waiting time. From these con¬ 
ditions it follows that H^g = H — ^ L*Li, where H is 
a selfadjoint operator interpreted as the system’s Hamil¬ 
tonian when isolated from the environment. 

After the first detection, the process is repeated start¬ 
ing with the state J[tw, ilf/WJitw, ijfW, with an indepen¬ 
dent waiting time. A full detection process is given by a 
finite measurement trajectory X = {{ti,ii), ..., (tmin)), 
where 0 < < • • • < . Each such trajectory has 

the final time T[X] = and total number of jumps 
A'[X] = n. When integrating over trajectories we will 
use the notation 


dX = 


OO 


E 


n—O ii,...,irL 


dti • • • dtj] 




( 2 ) 


Since we are dealing with quantum systems, it will turn 
out to be convenient to understand classical trajectories 
X as labels for quantum states |X) spanning an abstract 
output space of wavefunctions 


If ai{t) are the (singular) Bosonic annihilation operators 
at time t, and |H) is the vacuum state, then |X) = 
(^i) ■ ■ ■ conveys the in¬ 

tuitive idea that each detection event corresponds to a 
“particle” in the output. 

From the experimental point of view, there are two nat¬ 
ural ways of obtaining finite trajectories from the above 
scheme; the distinction has not been made explicit in 
the existing literature, and our first aim in the present 
paper is to clarify it. The first one repeats the single 
detection process K times. Clearly this is particularly 
relevant when the detectors have to be initialised after 
each detection, so that each measurement is actively per¬ 
formed on the system. This scheme has an associated 
state transformation given by, 

— 'db {tw.n Hn) ‘ J (tw,2 (f'u;,! : 

where = tn—tn-i, tw,i = ti are the waiting times be¬ 
tween the detection events, or quantum jumps. In words, 
given a trajectory X resulting from this process, the sys¬ 
tem is at the end in state J[X.]ip. 

The second scheme for trajectories is to let the process 
run until a given final time r. This suits better situations 
where the system is evolving in time, and detections take 
place spontaneously. The associated state transforma¬ 
tion in this case is, 

H”[X] = J[X]. (4) 

In this case, given a trajectory X resulting from the pro¬ 
cess terminated at time r, puts the system in the hnal 
state V'^['X.]ip due to the contractive evolution after the 
final jump. 

We now construct suitable output states describing 
these two processes, and relate them to the evolution 
of open quantum systems. 


A. Measurements with fixed number of jumps 

We fix a number of jumps K, and repeat the single de¬ 
tection process K times. For each realisation, we record 
the trajectory X, together with the corresponding condi¬ 
tional system evolution J[X]'0, and arrange them into a 
quantum superposition over all possible realisations: 


4- = y dX (Xl^-) |X), 

with orthogonality relations (X|X') = (5(X — X') match¬ 
ing with the integral ([^. Due to the ordering of the 
time points, the wavefunctions 4>(X) = (Xldt) corre¬ 
spond uniquely to symmetric functions of Ar[X] uncon¬ 
strained time variables on [0, oo). Hence, we may identify 
the output space with the Bosonic Fock space 

OO 

J-out 01,2 ([0,(X)),C^^)®'’" (3) 

n—0 


|MPSif)= [ dXJ[X]|V')(g)|X). (5) 

This pure state belongs to subspace of consisting 

of the tensor product between system and the AT-th layer 
of the output space ([^, and it is easy to check that it 
is normalised due to Its reduced density matrix on 
the output alone will be denoted by which is finitely 
correlated [28]. Its matrix elements are (X'|p|^®|X) = 
tr[pJ[X]’l'J[X']], where p = \'ip){tp\. For X = X', this is 
just the probability density for the trajectory X. 
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If we decouple the system at the end via suitable mea¬ 
surement (instead of tracing over it), we get a pure state, 
which is of matrix product form [29] with bond dimension 
at most the dimension of the system. This can be seen 
as follows: we first map the iiT-layer of the output space 
unitarily to the if-fold tensor product of the Hilbert 
space with waiting time basis *); where 0 < < oo, 

1 = 1 ,... Nl, by making a (unitary) identification 

|(^1:H)5 ■ ■ ■ 5 i^kY) — ; H) ^ ^ • 

We then choose a basis {|m)} for the system, and define 
normalizable output vectors 



Now clearly 

|MPSif) = ^ \mK){mo\tp}‘^ 

^niK ,mK-i ^ ‘ ■ ® ^m2 ,mi ^ ,?7io 5 (6) 

which can easily be put into an MPS form by choosing 
an orthonormal basis for the linear span of the vectors 
'I'm.m'- Note that this just amounts to diagonalising the 
associated Gram matrix 



We call ([^ the waiting time representation of the state 
|MPS^). 

It is also useful to find the unconditional evolution of 
the system under the measurement process, obtained by 
averaging over all realisations. This is simply obtained 
from |MPSif) by tracing over the output: 

tro,t[|MPS;f)(MPSK|]= f dXtr[J[X]pJ[X]t] 

JX:if[X]=if 

= Tf(p), 

where T* is the Schrodinger picture version of the quan¬ 
tum channel 

POO 

T(-)=^/ dt^J[t^,t]\-)J[tn,,i], 

describing one measurement step in the Heisenberg pic¬ 
ture. Since the process is discrete, the subsequent evolu¬ 
tion is obtained by iterating T. 

B. Measurements with fixed final time 

When the observation is terminated at fixed time r, 
each realisation consists of a trajectory X with T[X] < t. 


and the corresponding operator H’’[X]. We record this 
information into the state 

|cMPS^) := / dXH^[X]V'(8) |X). (7) 

./x:r[x]<T 

We recall that the output Fock space ® satisfies the 
factorisation property G where each 

term is a Fock space of trajectories up to and respectively 
after time r. Under this factorisation, the above state 
belongs to the subspace % ® j G II(t,oo) C 77 0 J"™* 
of system and output states with “spatial volume” re¬ 
stricted by the final time t. We note that the state can 
be written (see e.g. [ 20 ]) implicitly in terms of a path 
ordered exponential: 

|cMPS^) = 0 1 ^)), 

Conditional states obtained from this state by decou¬ 
pling the system via suitable measurements, are called 
Continuous MPS (cMPS), which were relatively recently 
introduced as useful variational states for quantum field 
theory [211132]. For later purposes we denote the output 
reduced state by p'f^. 

The averaged system evolution is now 

trout[|cMPS,)(cMPS,|] = f dXtr[I/”[X]pU”[X]t] 

7x:T[X]<t 
= e”W.(p), 

where the last equality is the Dyson expansion of the 
exponential with 

Nl 

wi-) = -ni-) + Y,Lli-)L,, 

and TZ{-) = i{-)HeS — iHlg{-). Of course, this is the Lind- 
blad generator of a continuous Markovian open system 
evolution (semigroup of channels) t i—>■ e'^'^. In physics 
literature, the measurement process is often called (quan¬ 
tum jump) unravelling of the evolution [22] . We have 
here presented it “the other way around” so as to em¬ 
phasise its similarities and differences to the above fixed 
K scenario, which develops more naturally from the jump 
picture, and has not been extensively investigated in the 
literature so far. 

From the mathematical point of view, the construction 
of the abstract output space is an instance of the mea¬ 
surement dilation theory of Naimark (cf. e.g. [21]) and 
Davies [22], which states that counting processes aris¬ 
ing from generalised (POVM) measurements can always 
be realised projectively on a larger Hilbert space. In 
concrete examples, the Lindblad generator often results 
from a unitary evolution coupling the system with some 
physical environment, typically a Bosonic heat bath. In 
this case, the abstract output field 7^°“* can naturally be 
identified with the environment, and the measurements 
themselves can be realised physically by measuring suit¬ 
able field operators. This leads to input-output models 
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for optical systems and stochastic Schrodinger equations 
[9] or filtering equations [U] describing the conditional 
evolution of the system monitored through the environ¬ 
ment. 


III. THERMODYNAMICS OF TRAJECTORIES 


Viewing the abstract output as a physical field with 
trajectories corresponding to configurations of excita¬ 
tions or “particles” on the output allows us to employ 
techniques from non-equilibrium statistical mechanics 
(see e.g. [35]) to study their properties at the limit of 
long observation times. The mathematical underpinning 
is the large deviations theory for the relevant extensive 
observables, which we now proceed to describe. 

We are interested in quantities obtained by increment¬ 
ing with some amount at the addition of each particle. 
Such a quantity is of the form F[X] = F{tw,n, *«), 

where F(tw,i) is a (possibly vector valued) quantity de¬ 
pending only on a single waiting time tw, and some prop¬ 
erty of the system we are monitoring, say “spin” i. Such 
quantity can be represented as the diagonal field operator 

F = J d^XF[X]|X)(X|, 


and the expectation values of any function f{F) in 
the above two processes have quantum expressions 
tr[p^f{F)] and tr[yo“^^/(F)]. In particular, by taking 
f{F) = S{F — Fq) gives the probability distributions of 
the observable F at point i^o- 

The most important time-extensive quantities are the 
particle number K and total volume T, obtained by tak¬ 
ing F{tu,,i) = 1 and F{tw,i) = respectively. These 
correspond to the above two ways of truncating an in¬ 
finite trajectory, leading to statistical ensembles of field 
particles with either fixed particle number K and fluctu¬ 
ating T, or fixed volume T and fluctuating K. In both 
cases we can also measure a “spin” operator M corre¬ 
sponding to F{tiu, i) = M{i), where M is some (in general 
vector valued) quantity depending on i. The associated 
probability distributions are given by 


Pk{T, M) := tr (^(T, M) - (T, M) 

Pr{K, M) := tr f {K, M) - {K, M) 


Our subsequent analysis rests crucially on the properties 
of the associated generating functions 

Zk{x,c) := tr[p|Pe-“^-“] = tr[pT^,(I)], 

Zr{s,c) := tr[pfCe-^^—*] = tr[pe”W-(I)], 

written in terms of the deformed generators [22] T^,,, and 
Ws,c, obtained from T and W by replacing J[X] and 
V”[X] with 


Jx.c[X] 


— Jx,c(j'W,n^ ^n) ‘ ‘ ‘ Jx^c(j^w,2-i '^2)Jx,c(j' W,1 •) ii), 


respectively, where Jx,c{tw,'i) = e 2 
Explicitly, we have 


T,,, = (xId-b7^)-l , 

Nl 

Ws,c = -n{A) + e-^ ^ (•)T^, 


( 8 ) 

(9) 


where the inverse 

poo 

(a;Id-f iR.)-i = / 

Jo 

exists whenever ||e“*d^fett-*a:/ 2 ) |j ^ which holds for 
all X > Xmin, where Xmin = 2ImAo < 0 and Aq is the 
eigenvalue of FleS with maximum imaginary part [43] . 
We restrict to a; > Xmin subsequently. 

We now suppose that has a unique eigenvalue 
gs(a:.c) equal iq Hg spectral radius, and that Ws_c has 
a unique eigenvalue 0{s,c) with largest real part. We 
further assume that both eigenvalues are nondegener¬ 
ate, associated to strictly positive eigenvectors F^^c and 
Es^ci respectively, and that they are the only eigenvalues 
with positive eigenvectors. We then assume that 
and 9{s,c) are also the dominant eigenvalues of the cor¬ 
responding Schrodinger picture generators (Ta;_c)* and 
(Ws,c)* in the same sense, with unique eigenvectors ax,c 
and CTs.c- These conditions hold under the generic as¬ 
sumption of strong irreducibility (primitivity) of the gen¬ 
erators |13|, but this is not a necessary condition (see 
Sec. 0. 

With the usual normalisation (see e.g. [55]) ir[ax,c] = 
tr[crs,c] = tr[Fx,cO-x,c\ = t.r[Es,cCrs,c\ = 1, we then have 


^-Kgix,c)^K^ tr[a,,c(-)]i"x.c (10) 

tr[a,.e(-)]l?..c. (H) 

In particular, the generating functions satisfy Zk{x, c) ^ 
gK g(x,c) Zr(s,c) ^ ^t these limits. The ex¬ 

tensive parameters A or r being large indicates thermo¬ 
dynamic limit at the output field. The above mathe¬ 
matical result that g(s, c) and 9(s, c) are well-defined in 
this limit lets us interpret them as thermodynamic poten¬ 
tials, and Zk{x,c) and Zx{s,c) represent path canonical 
partition functions [45], with the intensive variables x, 
s, and c acting as control parameters for the extensive 
quantities T, K, and M, respectively. In particular, we 
can identify x with pressure (of the output excitations), 
and s the chemical potential (controlling creation of new 
excitations). 

We proceed with the thermodynamic analogy by inter¬ 
preting the partition functions as normalisation factors 
for the corresponding path canonical ensembles, given by 
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FIG. 1: Sketch of trajectory ensembles. Time runs from left to right. Wiggly lines indicate quantum jumps, and straight ones 
deterministic evolution between jumps under Jfeff- (a) The s-ensemble. All trajectories in the ensemble are of a total fixed 
time r but can have any number of quantum jumps. Furthermore, between the time of the last jump and the final time there 
may be a period of no-jump evolution, (b) The a:-ensemble. All trajectories in the ensemble have the same number of quantum 
jumps K, but their time extension can fluctuate. In this case trajectories terminate after the iF-th jump. 


quantum states 


are 


:= Z;y(a;,c)-5e-5(^^+"'^)|MPSiy) 

= —L^[ dXJ,,,[X]|^)®|X) (12) 

Zk{x,c)2 J:x.-.k[x]=k 


|cMPS^“J := Z^(s,c)-5e-5 

1 f 


Zr(s,c)2 Jx:T[X]<t 


|cMPS^) 
dXP,^JX]V’® |X). 


(13) 


We denote the corresponding reduced output states by 
P^x c the corresponding classical trajectory 

ensembles (X|p^^^ j,|X) and (X|p(:^g®^|X) are sometimes 
called x-ensemble and s-ensemble, respectively [52]. 

From the asymptotic behaviour of the partition func¬ 
tions it follows by the Gartner-Ellis theorem that the cor¬ 
responding probability distributions satisfy the LD prin¬ 
ciple, when properly scaled: 


IMPS 


S{{T,M)-{T,M))\MPSk) 

K:T,M' ““ 


^Pk{T,M) 


'xgTk;- 


^Pk{T,M) 

j[XM(z\x) 


(16) 


where Tk-tm ■= {X : K[X\ = K,T[X\ = T,M[X] = 
M} is the space of trajectories with K jumps, the last 
jump occurring at time T, and the ’’spin” observable tak¬ 
ing value M, and 


|cMPS“5,,^) := 
1 

" '^Pr{K,M) 


5{{k,M)-{K,M))\cMPSr) 

kPr{K,M) 

f P”[X]V>®|X). (17) 


Pk{T,M) AT large (14) 

Pr{K,M) ^ r large. (15) 

These functions represent mierocanonical partition fune- 
tions for the above thermodynamical ensembles on the 
output field, determined by fixed extensive variables 
T, M, K. The associated thermodynamic potentials 
4>(t,m) and ip{k,m) only depend on the appropriately 
scaled intensive quantities t = TjK (total volume per ex¬ 
citation) and m = M jK (average value of M per excita¬ 
tion) for the ensemble of fixed K, as well ask = K/t (to¬ 
tal number of excitations per unit volume) and m = M/t 
(average value of M per unit volume) for the ensemble of 
fixed volume. The associated path mierocanonical states 


where Tr.KM ■= {X : r[X] < t,K[X\ = K,M[X\ = 
M}. Note that the mierocanonical states are singular 
vectors, when seen as functions over all trajectories, but 
can also be viewed as bona-fide Hilbert space functions 
over the restricted sets of trajectories, and are normalised 
in the latter sense. 

The potentials g and (j) are connected to each other by 
the Legendre transform, and so are 9 and (/?. Assuming 
that the functions are smooth and convex, we can apply 
the usual machinery of classical thermodynamics to get 

g{x, c) = —4>(t, m) — tx — m ■ c, (18) 

9{s, c) = —g}{k, rh) — ks — rh ■ c, (19) 
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where 

t{x, c) 
m{x, c) 
k{s, c) 
m(s, c) 



= lim 

{T[n)K,x 

dx 

>-00 

K 


= lim 

(M,[X])k.. 

dcj 

K—¥oo K. 

de 

= lim ■ 


ds 

r— >-oo 

T ’ 

do 

= lim 


dcj 

T —>-00 

r ’ 


s-ensemble 


fixed 

T = T*, any K, t* >tK 



K = K*, T* > tK- 


K = K* 

T=T* = tK- 



fixed K = K*, any t, t = tK- 


x-ensemble 


and the brackets denote expectation value with respect to 
the associated canonical ensembles. These familiar rela¬ 
tions describe the thermodynamics of the output field, in 
terms of the relevant intensive quantities obtained from 
the underlying continual measurement process. 


IV. EQUIVALENCE OF ENSEMBLES 

In classical statistical mechanics, the concept of equiv¬ 
alence of ensembles appears in different levels. Here we 
look at analogous statements for our quantum setup, the 
connection between the two different measurement sce¬ 
narios, fixed K and fixed r. 


A. Equivalence of potentials 


Often the equivalence is stated on the level of thermo¬ 
dynamic potentials. We have already seen that the Leg¬ 
endre transformations (181 and (191 establish one-to-one 


correspondence between the potentials for path canoni¬ 
cal and microcanonical partition functions, separately for 
systems of fixed K and fixed r. In [22], it was shown that 
there is also a one-to-one-correspondence between poten¬ 
tials g{x, c) and 6*(s, c). Since that result will be essential 
also for our subsequent development, we present it here 
in more detail. 

From (|^ and ([^ it follows that 


- end = e^xlA + n)-^ o {Ws,c - xId), (20) 


holding for all admissible values of the parameters. This 
connection between the generators is the fundamental 
reason for the ensemble equivalence. It implies that for a 
fixed c and a;, the eigenvector of Tx,c corresponding 
to the dominant eigenvalue is also an eigenvec¬ 

tor of Ws,c, with eigenvalue x, where s = g{x, c). By 
the uniqueness of positive eigenvectors (see Sec. jlllj ), it 
follows that 


X = 6{s, c), 

F ^ O' 


s = g{x,c), 
with some a. 


> 0 . 


( 21 ) 

( 22 ) 


Using the above thermodynamic relations, together with 
(211, we see that the intensive quantities defined via sys¬ 
tems of fixed A —)• 00 are related to those defined via 


FIG. 2: (Color online) Sketch of concentration of measure. 
The s-ensemble (in blue) corresponds to all trajectories with 
fixed overall time r* and any number of quantum jumps K, 
and where the time of the last jump tK is in general smaller 
than r. A subset of this set is that of trajectories where the 
number of jumps is exactly K*. The a:-ensemble (in red), 
in contrast, is that of trajectories all with fixed number of 
jumps K* but any time extent t where trajectories terminate 
with the last jump, r* = Ik* ■ The subset of trajectories 
with T exactly r* (in purple) also belongs to the s-ensemble. 
The two ensembles are equivalent when concentrated to this 
intersection if the fields s and x defining the ensembles are 
such that {K)s = K* and (t)^ = t*, see Sect. IV.B. 


systems of fixed r —)■ oo by 


t(x,c) = l/k{s,c), m{x,c) = rh{s,c)/k{s,c). (23) 

These relations show how the two different ensembles cor¬ 
respond to each other. In particular, the first relation in 
(23) states that T/K (total volume per excitation), for a 
given pressure x in the ensemble of fixed K, is equal to 
the inverse of KjT (number of excitations per total vol¬ 
ume), in the ensemble of fixed T, for a specific chemical 
potential s uniquely determined by x. 

By taking the adjoint of (20), we get the correspond¬ 


ing relation for the Schrddinger picture generators (T^, c) 
and (Ws,c)*- The dominant eigenvalues are again 
and 0{s,c), respectively (by our assumptions in Sec. Ill), 


so the same correspondence holds. However, the corre¬ 
sponding eigenvectors and Ps^c (normalised to having 
trace one) are not equal, but related via 


Ps,c = (rfd -I- A*) ^{px,c)/ trace . 


(24) 


B. Equivalence of concentration 

Another way of formulating equivalence of ensembles 
is in terms of the asymptotic behaviour of the ensem¬ 
bles (i.e. probability distributions) themselves. In equi¬ 
librium statistical mechanics, the equivalence between 
canonical and microcanonical ensembles states that by 
appropriately scaling the extensive variable (typically 
energy) held fixed in the microcanonical ensemble, the 
canonical ensemble will concentrate around that value in 
the thermodynamic limit. 






















In order to put our result into a context, we first 
briefly review the corresponding classical result, proved 
by Touchette [20], using the above “thermodynamics 
of trajectories” approach. (For similar results, see e.g. 
[201 [13 mi)- Given an observation time r, let X,- de¬ 
note the corresponding time truncation of any infinite 
trajectory X, and suppose we have a family of path 
probabilities Pt-(Xt-), r > 0 specified by measurement 
process on some nonequilibrium classical system. Con¬ 
sider the path microcanonical ensemble P™5f(XT) := 
5{K — K[Kr\)^TO^T)/P t{K), where is Pt{K) assumed 
to satisfy the LDP, i.e. Pt{K) '■= f d'KrS{K[X.r] — 
iF)PT-(XT-) ~ Then the path canonical ensem¬ 

ble P“"(X^) ;= e-®-^[^-lP^(X^)/Z^(s) has LDP with 


Zt-{s) := 6 “^PriK) ~ where (19) connects 


the potentials, we set k{s) := —dO/ds, and Kr := rfc(s) 
for fixed s. Then 


|(MPS5,“^J^((r,M) - iTK,MK))\cMPS 


can \ 1 2 
tk\s,c) \ 


Zk{x,c) ■ Zrj^(s,c) 




K large, (26) 


This shows that up to sub-exponential corrections, both 
path canonical states concentrate on the common “mi¬ 
crocanonical” shell where the extensive variables T, K, 
M have their typical values determined by the intensive 
quantities t{x,c), k{s,c), and m(x,c). Note that there 
are only two free parameters here: we can fix either x 
(pressure) and c (in which case s = g{x,c)), or we fix s 
(chemical potential) and c (in which case x = 9{s,c)). 
The relations (23) ensure that the limits iL —oo for sys¬ 
tems of fixed K, and r —>■ oo for systems of fixed r are 
consistent. 


P“"(X,)/P“5, (X.) = e° 


for large r. 


(25) 


C. Operational equivalence 


for those infinite trajectories X whose time truncations 
X,- satisfy KlK-r] = Kr for each r. Here o(r) denotes 
sublinear growth in r, i.e. lim,- o{t)/t = 0. In words, this 
result states that the canonical ensemble “concentrates” 
on the microcanonical shell with iF[X,-] = rfc(s) [IH] - 


In our quantum case, we can formulate an analogous 
statement in terms of the scalar products of the rele¬ 
vant quantum states. Hence, using the definitions (16) 
and (0, the equivalence between the path canonical and 
microcanonical ensembles reads 


|(MPS^^“ 

|(cMPS“"JcMPS“5,^,m.)I" = 


The above concentration equivalence only considers 
the overlap of the state on the common singular support 
specified by the typical values of the intensive quanti¬ 
ties. We now proceed to formulate an operational equiv¬ 
alence stating that full statistics of all counting measure¬ 
ments up to a fixed time tq, are the same for both states 
and |cMPS):“,,), in the thermodynamic limit 
iL —>■ oo and r —>■ oo. 

More precisely, we consider the scenario of long time 
(r —>■ oo) and large number of jumps {K —>■ oo) evo¬ 
lution, but only observe the output state up to a fixed 
time To and compare the reduced output states in the 
two ensembles. The restriction correspond to the tensor 
product factorisation of the output Fock space 


where Tk = t{x,c)K, Mk = m{x,c)K, Kr = 
ko{s,c)T, and Mr = rh{s,c)T. Both results can be 
proved easily from the definitions, using the above es¬ 
tablished LDPs. As noted before, the micro canonical 
states K[^) are singular vectors in the space 

of all trajectories with fixed K, but the inner product 
can be computed by interpreting the state as a delta dis¬ 
tribution. Similarly, we can prove that the equivalence 
of the type ( p^ ) holds for the classical path probabili¬ 
ties Fk{Xk) := (X|p|P|X) and P,(XO := {X\pf^\X) 
corresponding to the fixed observables T and K, respec¬ 
tively. 

Building on the nontrivial connection (21), we now 
investigate a form of asymptotic equivalence between 
the canonical discrete and continuous time ensembles 
and respectively |cMPS):“j,). For given pairs 
(x, c) and (s,c) related as in mi)’/ we define tk '■= 
Kt{x,c) = K/k{s,c) and Mk ■= Km{x,c), where K 
is considered to be large. If r := tk, then the two en¬ 
sembles have the same means of extensive observables. 
Moreover, the quantum ensembles can be compared when 
restricting to trajectories for which T(X) = tk, i.e. the 
evolution P” after the last jump is trivial, cf. (|^. 


= PoZ)P+, (27) 

where the left (right) terms are Fock spaces of the count¬ 
ing fields over (0,to] and (to,oo), respectively. This cor¬ 
responds to the identification |X) = |Xo) G |X+) of basis 
vectors, where any finite trajectory is split into concate¬ 
nation X = Xq V X+ so that Xq is the part with all 
jump times at most tq, and X+ the rest. The trajectory 
integral separates as / dX = f dXg f dX^. 

Given an arbitrary state g on the full output its 

reduced density matrix to Pq reads 

(Xo|tr+[^.]|X') = y'dX+(XoVX+|e|X' VX+), (28) 

where tr+ is the partial trace over JF|_. We will also need 
the projection onto the A^th ’’layer” of the Fock space Pq 
which is given by 

P^= f dXo|Xo)(Xo|; 

J{Xo:K[Xo]=N} 

these are just the eigenprojections of the number opera¬ 
tor X]w=o ^ defines the basic counting measure¬ 

ment for the time interval [0, rg]. 
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We need to determine the reduced states tr_|_ ^ 

and tr_|_ . Note that for s = c = 0 we have 

tr+[p^^g^g] = i.e. the state describing the physi¬ 

cal evolution on [0,rg] remains the same for all t > rg, 
as it should due to the temporal causality of the physical 
setup. However, the biased state tr_|_ does not have 

this property^ as we will see explicitly from the compu¬ 
tation below. (The same phenomenon happens also in 
the classical case, as pointed out e.g. in [36].) The re¬ 
duced state tr+ J is of course not expected to be 
remain constant with increasing K even for the physical 
X = 0 case, because in this case the jump times are not 
temporally constrained. 

We will now proceed to prove the following result which 
holds for given (s, c) and x = 9(s,c): 


explicitly how the Dyson expansion of the deformed “evo¬ 
lution” factorises along this split. For the last part, 


g(r-ro)W„. ^ J dx+F+/[X+]t(.)Ky[X+]. (32) 

The matrix elements of the total Fock state are 


(xip^pc |X') = ^ 


tr[pF-JXfF-JX]]^ 


Zris,c) 


hence using (28) and (321, we get the reduced state: 
(Xo|tr+[p^fi;J|X') 

_ tr[pF”oJX'o]te(”-”°)w^.-(I)F”g[Xo]] 
Zr{s,c) 




T —>-00 oo 

-^ pZ- 


K—¥oo 




iV=0 




,oo pu 
s,c-^ 


(29) showing that 


(30) 


tr+[p 


cFC 1 
r;s,cJ 




g(T-ro)(W,,c-e(s,c))^][^ 

gToe(s,c)^-Te{s,c)zZs, c) 


S.. 


where the convergence is in the trace norm of the outpuf 


Ps.c ^ 


— p-ToS(e,c) 


and S', 




(O.ro) 


tvlpEs^c] ^Sl^Es^cSs,c, (31) 

H is the operator defined by 
Ss^c|Xg) = Fj,'(°[Xg]'(/;, for all Xg. Note that p^^ in (29) is 
a finite-rank operator that lives in the space spanned by 
the vectors S| ^m), where |m) are the basis vectors of the 
system Hilbert space. In contrast, the limit state in (30) 
has infinite rank because the projected states 
are all nonzero. 


The interpretation of (29) and (30) is that in the long 


run the reduced states of the two ensembles converge, and 
the limits have the same diagonal blocks with respect to 
the decomposition of the Fock space into fixed-lV lay¬ 
ers. From the operational point of view, this means that 
in the limit, the two ensembles yield the same measure¬ 
ment statistics for any quantum observable on Eq that 
commutes with the number operator 
compatible with the counting measurement). However, 
while the time ensemble state may exhibit coherences be¬ 
tween different layers, the discrete time ensemble state is 
block-diagonal. This inequivalence on the level of off- 
diagonal blocks is discussed in more detail in the next 
section. 

We start the proof by noting that the second limit is 
crucially based on the equivalence of potentials (211, and 
the eigenvectors (22). We first define 


F+/[X+] = Jj;),[X+], 


where is just with the first jump time ti re¬ 

placed by ti — Tg so as to shift the starting point to 
Tg where the trajectories X+ begin. (In case of an 
empty trajectory, we omit Js^c and put r[X] = 0 and 
T[X+] = Tg). According to the split ( |^ , we have 
F”^[X] = F+j”[X+]F”«JXo], for any t > tq. This shows 


Here Tg is constant, and the only r-dependence is in the 
system operator sandwiched in the middle. Since that 
operator converges to pg°(s,c) due to 0. and the sys¬ 
tem is finite-dimensional, we get ( |2^ . Note that for the 
case s = c = 0, we just have (I) = I, so the 

unobserved part of the trajectory does not contribute at 
all (as expected). For the biased case, we have to take 
the limit to get a fixed density matrix. 

We now proceed to prove the second limit (|30|). We 


fix a pair Xg, Xg which determines an element of 
the reduced density matrix of the full state p^^ 

If max(Ar[Xg], ATjXg]) < K, the relevant matrix ele¬ 
ments of the total output state are of the form (Xg V 
^+IPk%.cI^o V X+). Since p^^^^ lives in the ATth 
layer of the total Fock space, the element is zero un¬ 
less if[Xg] -b K[X+] = ATjX'g] -b X[X+] = K, i.e. 
ATjXg] = ATjX'g]. Similarly, if max(X[Xg], A:[X(,]) > K 
then the matrix element is zero since the full state has ex¬ 
actly K excitations. If max(Ar[Xg], ATjXg]) = K then all 
counts happen before rg and the contribution to the re¬ 
duced state is zero unless AljXg] = A'[Xq] = K. Putting 
all the cases together we see that the reduced state is 
block-diagonal with respect to the decomposition of Eq 
in Fock-layers. 

We will now fix = ATjXg] = Ar[XQ] and show that 
the Wblock of the reduced state converges to the cor¬ 
responding block of the limit state state (30). 

Since AT —>■ oo we can assume that N < K so the nonzero 
elements have T[X] = T[X+], and are given by 


(XgVX+|p|.%,JX'gVX+) 

tr[pFg;°[X;)]t(J”;^[X+]tj-;),[X+])Fg-»[Xg]] 
e^^°Z k{x, 


Note that the part depending on Xg and Xg is again of 
the same form as before, except that s does not appear. 
We can fix that by defining s = g{x,c), and noting that 

















10 


Vy^“[Xo] = [Xq]. We then use (28) to integrate 

over X+; now the integration only goes over trajectories 
with it'[X+] = K—N. By making the change of variables 
ti I—>■ ti + To on X+ we can transform the integral: 


In particular, each transformation only changes one of 
the ensembles. Note that the imaginary shift in the s- 
parameter changes the matrix elements of accord¬ 
ing to 




/{X+:iC[X+]=if-Ar} 


^x(T[y.+ ]-ro) 


J{X-.KIX\=K-N} 


dX J[X]^ 


and hence 


P'^tr+[pZ^x,c\PN — 




Now we use the fact that s = g{x, c), together with (101, 
to take the limit df —>■ oo; the system operator in the 
middle converges, and we get 


r„FC 1 dO if->00, oO ^l,c^x,c^s,c ^0 


P ^ tr + [pZ!^x,c]^n 


again in the trace norm. Since x = 0{s,c) and = 
Cis,cEs,c by (22), the operator in the middle coincides 
with p3°(s,c). To obtain (30) we note that the separate 
convergence of the Wblocks implies the convergence in 
trace norm of the full state, which can be shown by using 
a quantum version of Scheffe lemma, cf. e.g. Theorem 
2.20 in [37]. 


D. Inequivalent aspects 


Having established that the above three levels of “en¬ 
semble equivalence” hold, we close this section by em¬ 
phasising that the two ensembles are not equivalent in all 
aspects. In fact, from the statement of operational equiv¬ 
alence, we can already observe that there is be a differ¬ 
ence in coherences between different Fock layers; while 
the continuous time ensemble contains coherent super¬ 
positions of output states with different number of ex¬ 
citations, the discrete time output is a mixture of such 
states. A consequence of this is that the two ensembles 
can exhibit different properties, when the coherences are 
relevant. In order to see that this is indeed the case, 
we look at the following simple transformations of the 
process parameters H and Lp. 

(PI) H ^ H, Li^ e^'^Li. 

(P2) H ^ H + Li^ U. 

Here (/) is an arbitrary real number. It is straightforward 
to verify that these correspond to the following transfor¬ 
mations of the cc-and s-ensemble states: 


PI) pZx,c invariant, pfZ ^ 

P2) pfZc invariant, p^^x.c ^ pZi,x- 2 i 4 ,,c- 


(X|p“ ^JX') ^ e 


K'-J 




hence introducing the coherences not present in the diag¬ 
onal elements; hence this inequivalence is in accordance 
with the above equivalence results. The inequivalence of 
the ensembles in phase transformations (PI) and (P2) 
has consequences for parameter estimation of Markov 
processes; this is however beyond the scope of the present 
paper, and will be analysed in a separate publication. 


V. RENEWAL PROCESSES 

As already mentioned, the discrete MPS associated to 
the a;-ensemble is structurally simpler than the continu¬ 
ous MPS describing the s-ensemble. In order to illustrate 
the difference, we now consider a physically relevant ex¬ 
treme case where the former is a product state. 

The measurement scheme described above is a re¬ 
newal process, if there exists a specific state |0) such 
that Li = |0)(:/3i| for some vectors ipi. This means after 
each measurement, the system is in state |0), regardless 
of the measurement outcome. The jump operators are 

and the waiting time repre¬ 
sentation of the MPS is 

|MPSk) = |0) (g) <g) • • • (g) <g) 


where 


Nl 

fko,m = 'y 

i=l 

and |m) is a basis of the system Hilbert space con¬ 
taining the special vector |0). Hence, disregarding the 
first factor related to the initial state of the system, we 
simply have a tensor product of the single wavefunc- 
tion 5'o,o- The transition generator is trivial, T(-) = 
I(0|(-)|0), and the deformed generator is likewise not 
ergodic: T 2 ,(-) = (xld + 7?.)“^(II)(0|(-)|0), where D = 
— ~2Imi7eff. However, T^, clearly has only 
one nonzero eigenvalue = (0|(a;ld -I- 7l)~^{D)\0). 

Hence, the dominant eigenvectors px = |0)(0|, and 
Fx = (ccid -|- Tl)~^{D)e~^^^^ are uniquely determined, 
and we see directly that the limit ( [l0| holds. Moreover, 
p, = (a:Id + 7^*)-l(|0)(0|)/tr[(rfd + 7^*)■^(|0)(0|)]. 

A basic example of a renewal process is the sim¬ 
ple 3-level system with the ground state |0), and two 
exited states |1) and |2), with the transitions |0)-|1) 
and |0)-|2) driven by two resonant lasers with Rabi fre¬ 
quencies Hi and H 2 , respectively. In addition, level 
|1) decays to |0) with rate k, and we detect the pho¬ 
tons emitted in this process. Hence the Hamiltonian is 


dtyj IpPi 




to) \txu,i), 
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H = X]m=i ^m(|0)(w| + |w)(0|)j and we only have one 
Lindblad operator L = •yK|0)(l|, i.e. = 1- One 
can easily check that i7eff has three distinct eigenvalues 
with strictly negative real parts, except for some spe¬ 
cial values of the parameters; hence the normalisation 
conditions Q are satisfied with the single jump oper¬ 
ator J[tw) = •y/K|0)(l|e“*‘^'’®, which maps everything 
to the ground state. This example was used in m to 
demonstrate how intermittency in the jump trajectories 
can be related to dynamical phase transitions in the sys¬ 
tem. From the point of view of the present paper, the 
example is interesting because the two MPS states intro¬ 
duced above exhibit very different aspects of the dynam¬ 
ics. 


VI. CONCLUSIONS 

The “thermodynamics of quantum trajectories” for¬ 
malism m has been successfully employed in uncover¬ 
ing and analysing dynamical phase transitions in open 
quantum systems, through the statistical properties of 
an appropriately biased “s-ensemble” of quantum jumps 
trajectories. In [21] it was shown that the large devia¬ 
tions rates of the fixed-time s-ensemble can be deduced 
from those of an alternative “cc-ensemble” containing bi¬ 
ased trajectories with fixed number of jumps but random 
time length. 

Here we have strengthened this correspondence by lift¬ 
ing the ensemble equivalence to the level of the quan¬ 
tum states m of the two corresponding input-output 
systems: a continuous-time one governed by a deformed 
Lindblad generator, and a discrete one whose transition 
operator encodes the jump-to-jump dynamics. As a first 
equivalence result we showed that the overlap of the con¬ 
tinuous and respectively corresponding discrete system- 
output (matrix product) states decreases subexpontially 


in the limit of large times and counts. Furthermore, the 
restrictions of the output states to a finite time interval 
converge, and the limits have identical diagonal blocks 
with respect to decomposition of the Fock space in “lay¬ 
ers” of fixed number of excitations. However, a specihc 
quantum feature which is not present in the classical set¬ 
up [ini [H] is that the s-ensemble may exhibit coherences 
between different layers, while the a:-ensemble is diago¬ 
nal. A consequence of this inequivalence is that certain 
dynamical parameters may be identified by output mea¬ 
surements in one ensemble but not in the other. For 
instance, the off-diagonal blocks of the fixed-time output 
contain information about the phase of the jump opera¬ 
tor, while the fixed-jumps state does not depend on it. 

On the trajectories level, in certain situation the x- 
ensemble |22j can be more amenable than the s-ensemble 
[19] to numerical simulations via paths sampling tech¬ 
niques [35] • In the same spirit, we speculate that the 
above quantum equivalence results can contribute to¬ 
wards a better understanding of dynamical phase tran¬ 
sitions on the level of quantum states. A possible appli¬ 
cation is the extension to continuous time of the Sanov 
Theorem for the empirical measure of multiple successive 
jumps, developed in [119] . In a different direction, the two 
ensembles set-up could be used to unify the existing sys¬ 
tem identihcation and asymptotic normality theory for 
discrete [20] and continuous jH] quantum Markov pro¬ 
cesses. 
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